Halo Clustering with Non-Local Non-Gaussianity 
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We show how the peak-background split can be generalized to predict the effect of non-local 
primordial non-Gaussianity on the clustering of halos. Our approach is applicable to arbitrary 
primordial bispectra. We show that the scale-dependence of halo clustering predicted in the peak- 
background split (PBS) agrees with that of the local-biasing model on large scales. On smaller 
scales, k > 0.01 h Mpc -1 , the predictions diverge, a consequence of the assumption of separation of 
scales in the peak-background split. Even on large scales, PBS and local biasing do not generally 
agree on the amplitude of the effect outside of the high-peak limit. The scale dependence of the 
biasing — the effect that provides strong constraints to the local-model bispectrum — is far weaker for 
the equilateral and self-ordering-scalar-field models of non-Gaussianity. The bias scale dependence 
for the orthogonal and folded models is weaker than in the local model (~ A; -1 ), but likely still 
strong enough to be constraining. We show that departures from scale-invariance of the primordial 
power spectrum may lead to order-unity corrections, relative to predictions made assuming scale- 
invariance — to the non-Gaussian bias in some of these non-local models for non-Gaussianity. An 
Appendix shows that a non-local model can produce the local-model bispectrum, a mathematical 
curiosity we uncovered in the course of this investigation. 
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I. INTRODUCTION 

An active effort to seek departures from Gaussianity in 
primordial cosmological perturbations is currently under 
way [1]. The motivation stems from the desire to learn 
more about inflation [2, 3]. The simplest single-field slow- 
roll (SFSR) inflation models predict that any departures 
from Gaussianity should be undetcctably small [4] . How- 
ever, theorists generally consider these to be toy models, 
or effective theories, and there should be, at some point, 
phenomena that depart from the predictions of these sim- 
plest models. Many physics models for inflation [5-7], or 
alternatives/additions to inflation [8, 9], do in fact make 
predictions for non-Gaussianity at a much higher level 
than those of SFSR inflation. 

It has been proposed to search for primordial non- 
Gaussianity through a variety of techniqes, including 
cluster, void, and high-rcdshift-galaxy abundances and 
properties [10, 11]. Primordial non-Gaussianity had long 
been sought with large-scale-structure surveys [12], but 
the sensitivity of these surveys had generally been weak 
compared with that due to the CMB [13-15]. It has 
recently been shown, however, that the effects of pri- 
mordial non-Gaussianity may be amplified in the clus- 
tering of dark-matter halos [16, 17, 19, 20]. The number 
density of halos is determined by the threshold overden- 
sity for gravitational collapse. In a region with a long- 
wavelength overdensity, the upward fluctuation required 
for collapse is reduced and the abundance of objects 
thus higher. However, in non-Gaussian models, the local 
power spectrum (or density-fluctuation amplitude) may 
also be affected by the local long- wavelength overdensity. 
If so, this provides an additional dependence of the lo- 
cal abundance of halos on the long-wavelength density 



field. In the local-model of non-Gaussianity, this leads to 
a rapid increase in the biasing of halos on large scales. 
Null searches for this large-scale enhancement place con- 
straints to primordial non-Gaussianity that are compet- 
itive with those from the CMB [17, 21]. 

Still, the local model is just one of many possible mod- 
els for non-Gaussianity, and it is natural to inquire how 
other forms of non-Gaussianity, described by other bis- 
pectra, may affect biasing. This question was tackled 
in Rcf. [22] using local Lagrangian biasing (LLB) [23] to 
model the clustering of halos. However, different analytic 
descriptions of biasing that make similar predictions at 
lowest order may differ at higher order [24], and so it 
is important to investigate different biasing models to 
understand better the model predictions and theoretical 
uncertainties. In this paper, we calculate the clustering 
of halos for a general primordial bispectrum using the 
peak-background split (PBS). Our approach generalizes 
that in Rcfs. [16, 17] (see also [18]) and complements the 
local-bias calculation of Ref. [22]. 

As examples, we apply our results to the equilateral 
[25], folded [26], orthogonal [27], and self-ordering-scalar- 
field (SOSF) [9] bispectra, as well as the local-model bis- 
pectrum [13, 28]. We find that the scale-dependence of 
the biasing — the effect that leads to strong constraints 
to the local-model bispectrum where the bias correction 
goes as k~ 2 — disappears for the equilateral and self- 
ordering-scalar-field models. While the scale dependence 
for the orthogonal and folded models is weaker than in 
the local model (~ fc _1 ), it may still be strong enough 
to be constraining. We also find that small departures 
from scale-invariance in the primordial power spectrum 
may lead to large corrections to the non-Gaussian bias 
for some of these non-Gaussian models, especially for the 
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equilateral and SOSF models. 

We show that while the PBS prediction for the scale- 
dependence of halo clustering for general bispectra agrees 
with the LLB predictions at large scales, they differ at 
small scales, a breakdown of the separation of scales as- 
sumed in the PBS. Furthermore, the amplitude of the 
non-Gaussian halo bias is determined by different quan- 
tities in both approaches, and these quantities do not in 
general agree outside a specific high-peak limit of local 
biasing. 

The paper is structured as follows: we review the case 
of general, quadratic primordial non-Gaussianity in § II. 
§ III presents the peak-background split argument, and 
§ IV derives the non-Gaussian halo bias in this approach. 
§ V contains some results on the squeezed limit of differ- 
ent primordial bispectra. Finally, we compare the PBS 
and local-biasing approaches in § VI, and conclude in 
§ VII. An Appendix shows that the local-model bispec- 
trum can arise from a non-local model, a mathematical 
curiosity that arose during the development of the for- 
malism in this paper. 



II. GENERAL QUADRATIC 
NON-GAUSSIANITY 



A. Basic formalism 

Suppose the primordial potential $ is a general, non- 
local quadratic function of a Gaussian field 4>. Expressed 
in real space, the most general expression preserving sta- 
tistical homogeneity and isotropy is a two-dimensional 
convolution, 

$(x) = 0(x) + / NL j d 3 y y>ziy(y,z)0(x+y)0(x + z), 

(1) 

where the kernel W(y, z) is only a function of y, z, and 
y • z and is symmetric, W(y, z) = W(z,y). Here, the 
non-Gaussianity amplitude is quantified by the parame- 
ter /nl- Although the spatial average of $ is not zero, 
but of order {(f> 2 ), this constant offset is unobservable and 
can be eliminated by a gauge transformation. 
In Fourier space, Eq. (1) reads 



$(k) = 0(k)+/ NL 



d 3 fci 



(2tt) 3 

x Wpcx.k-kxMkiMk-kx), 



(2) 



where W is the Fourier transform, 



d 3 y / d 3 ze- iki y - ik2 Z M/(y,z). (3) 



W(ki,k 2 ) 



Statistical homogeneity requires that W has no depen- 
dence on the directions of ki and k 2 , only on their relative 
directions. It is thus a function only of the magnitudes k\ 
and fc 2 (and symmetric in these arguments) and the dot 



product ki • k2. It can thus alternatively be written as 
a function W(ki, k 2 , fc 3 ) of the magnitude fc 3 = |ki + k 2 | 
of the third side of the triangle constructed from kx and 
k 2 . 

The leading non-Gaussian correction to the potential 
is commonly parametrized in terms of the bispectrum 
B(ki, fc 2 , fc 3 ) of $, defined by, 

($(kx)$(k 2 )$(k 3 )) = (27r) 3 ,5z 5 (k 1 +k 2 +k 3 )S $ (fc 1 ,fc 2 ,fc 3 ), 

(4) 

where 5d is a Dirac delta function. It is straightforward 
to calculate the bispectrum of $ given Eq. (2). We obtain 

fca, As) = 2/nl| W ? (kx,k 2 )P*(A;i)P 4 (fc 2 ) 



-2 perm. 



(5) 



Here, P$(fc) denotes the power spectrum of <I>, which to 
leading order in /nl agrees with that of <f>. The two 
permutations not written are the two remaining cyclic 
permutations of k\, fc 2 , fc 3 . The requirement of statistical 
isotropy dictates that P$ is a function only of the trian- 
gle shape, not its orientation, and so the bispectrum can 
be written as a function B$(ki, fc 2 , fc 3 ) only of the mag- 
nitudes of its three arguments. The Fourier-space kernel 
W is, on the other hand, only required to be symmetric 
under exchange of its two (vectorial) arguments. Hence, 
Eq. (5) does not uniquely define W; there may be sev- 
eral different W that yield the same bispectrum (note 
however that these will give different trispectra at order 

/nl)- _ 

The ambiguity in W can be eliminated if we assume, 
for example, that it is symmetric under exchange of its 
three arguments. With this additional assumption, the 
Fourier-space kernel can be written in terms of the bis- 
pectrum as 



Wsym(fcx,A; 2 ,fc 3 ) 



1 



B^(ki,k 2 ,k 3 ) 



2/NLPP2 + P1P3 + P2P3 



(6) 



where Pj = Pj>(fci). Eq. (6) becomes more clear by con- 
sidering concrete examples. In the local model, the bis- 
pectrum is given by [4, 13, 28] 

B% c (k 1 ,k 2 , fc 3 ) = 2./nl [P$(fci)P$(fc 2 ) + 2 perm.] . (7) 



Applying Eq. (6) immediately yields 

and going back to real-space, we obtain 
W loc (y, Z ) =6 D (y)6 D (z). 



(8) 



(9) 



This, together with Eq. (1), yields the well-known local 
expression for $, 



$(x) =0(x) + / NL </> 2 (x). 



(10) 
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Finally, we note that the prediction of the scale- 
dependent non- Gaussian halo bias on large scales does 
not depend on the particular choice of kernel. As we will 
see in § IV, this prediction depends on the squeezed limit 
of the kernel, where fc 2 3> k\. Denoting q = fci/fc 2 , this 
limit corresponds to q — > 0, and in the same limit Eq. (5) 
can be written as 1 



w(ki,k 2 ; 



squeezed 



>(fci, k 2 , k 3 ) 



4/N L P$(fci)P*(fe) 



(1 + 0(9)) 



fill 



Thus, in the squeezed limit the kernel W and in par- 
ticular its leading scaling with q arc uniquely specified. 
Furthermore, Eq. (5) also shows that all kernels agree for 
equilateral triangles, further restricting the magnitude of 
the kernel-choice ambiguity. In § V we will return to 
the squeezed limit of W and show how it determines the 
scale-dependence of the non-Gaussian halo bias. 



B. Some primordial bispectra 

We now summarize other non-local bispectra that have 
been discussed in the recent literature and that we will 
consider here. All of these bispectra are scale-free [i.e. 
they scale as Pj under a uniform rescaling of all ki] , but 
we reiterate that this is not necessary for Eqs. (2) and 
(6) to apply. These bispectra simplify considerably in the 
squeezed limit, which we will discuss in § V. 

In the equilateral model [25], the bispectrum is given 
by 

Bf = 6/ NL [(-PiP 2 - 2 perm.) - 2(P!P 2 P 3 ) 2 / 3 

+ (Pl /3 P^ /3 P 3 + 5 perm.)] . (12) 
The folded bispectrum is given by [26] 

Pi 1 = 6/nl [(P1P2 + 2 perm.) + 3(P!P 2 P 3 ) 2 / 3 

_(pV3p2/3p 3 + 5perm) ] (13) 

The orthogonal bispectrum is given by [27] 

i/ NL { - 3 (PiP 2 + 2 perms) - 8(PiP 2 P 3 ) 2/3 



P£' th = 6j 



(PiP 2 1/3 P 3 2/3 ) + 5pcrms }. (14) 



The bispectrum due to self-ordering scalar fields (SOSFs) 
is [9] 



R SOSF 



437 / NL (P1P2P3) 
k 3 



2/3 



1 



0.485^ 
fci 



1 



1.87-^ 

fci 



( 7^ 

M 



(15) 



1 This assumes that the kernel in the last permutation of Eq. (5), 
involving the two large k values, does not grow faster than <?~ 3 . 



where in this case we restrict k\ > fc 2 > k 3 . We adopt the 
usual convention of defining /nl so that P$ = 6/nl(P*) 2 
for equilateral triangles (fci = fc 2 = k 3 ), except for the 
folded bispectrum, which is zero for equilateral triangles. 



C. Relating the processed and primordial power 
spectra 

So far we have been discussing the primordial poten- 
tial and its power spectrum and bispectrum, but what we 
need for galaxy clustering is the processed potential. The 
processed potential, which we denote with a subscript 
"0," is related to the primordial potential in Fourier space 
via $o(k) = T(fc) ( E > (k), where T{k) is the transfer func- 
tion. If we similarly define a processed Gaussian random 
field </>o(k) = T(fc)</>(k), then the processed potential may 
be written, 



$ (k) = <Mk) + / NL 



d 3 h 



(2tt) 3 

x W ; (k 1 ,k-k 1 )^ (k 1 ) ( /) (k-k 1 



(16) 



in terms of a processed kernel, 



wyki,k 2 ) 



r(|ki + k 2 |)ty(ki,k 2 ) 

T(k!)T(k 2 ) 



(17) 



Note that even if we choose W(ki, k 2 ) to be a symmet- 
ric function Ws ym (fci, fc 2 , k 3 ) of the three magnitudes fci, 

fc 2 , k 3 , the processed kernel, Wo(fci,fc 2 ,fc3) is not sym- 
metric under exchange of all three arguments. Still, it 
is symmetric under exchange of its two vectorial argu- 
ments if written Wo(ki,k 2 ). The configuration-space 
kernel Wo(y, z) is then obtained from the Fourier trans- 
form of W / o(ki,k 2 ). It is this processed kernel that will 
appear in the analysis below. 

Alternatively, one can define the processed kernel Wq 
directly from the bispectrum and power spectrum of the 
processed potential $0: m analogy with Eq. (6): 



Wo(ki,k 2 ) 



1 flg (fci,fc 2 ,|ki + k 2 |) 
2/nl P*o (fci )P* (fe) + eye. 



(18) 



where 



P* (fci, fc 2 , fc 3 ) = T(k 1 )T(k 2 )T{k 3 )B^(k 1 , fc 2 , k 3 ), (19) 
P* (fc) =T 2 (fc)P$(fc). (20) 

Note, however, that if we choose the processed kernel to 
be symmetric in all three arguments, then the primordial 
kernel will no longer be symmetric. That is, either the 
primordial or the processed kernel may be fully symmet- 
ric, but not both. 

In our analysis below, we will consider these two forms 
for Wo, one (which we label with the subscript "1") in 
which the primordial kernel is fully symmetric, and one 
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(which wc label with the subscript "2" ) in which the pro- 
cessed kernel is symmetric. To be clear, the two kernels 
we consider below arc 



Woi(ki,k 2 ) 



niki 



T(h)T(k 2 



sym(kl, k 2 ) 



and 



Wo 2 (ki,k 2 ) 



1 flgofo.fcajki+kal) 
2/nl P® (h)P$ (k2) + eye. 



(21) 



(22) 



There are, of course, other kernels one could choose. 
However, with any definition, Wq = W whenever all 
three k vectors are sufficiently small so that T(k) ~ 1. 
Moreover, they have the same squeezed limits: 



Wo(ki,k 2 ; 



1 



W(ki,k 2 ) 

g$(fci,fc 2 ,fc 3 ) 
T(fci)2/ NL P$(fc 1 )P$(fc 2 )' 



1 



(23) 



where we have used Eq. (11) for the last equality. This 
shows that the processed kernel is uniquely specified in 
the squeezed limit as well. Further, W i and W 02 also 
have the same equilateral limits, which further restricts 
the numerical uncertainty associated with this theoret- 
ical ambiguity. Finally, note that Wo ^ 1 in general 
for the local model, even though W is unity. This re- 
flects the fact that the relation between the primordial 
and processed potentials (and matter perturbations) is 
non-local. 

Before proceeding, we note that one consequence of our 
approach is the realization that a non-local inflationary 
model can give rise to the local-model primordial bispec- 
trum. This is spelled out in the Appendix. 



III. PEAK-BACKGROUND SPLIT FOR 
GENERAL QUADRATIC NON-GAUSSIANITY 

The physical, non-Gaussian, matter-density perturba- 
tion 5 is related to the processed potential $o through 
the Poisson equation, 



V 2 

— <f> = AttG pS. 



(24) 



Similarly, we can define a fictitious Gaussian density field 
5 related to the Gaussian potential 0o by the same rela- 
tion. In Fourier space, these relations read, 



3 (1 + z)H%VL 

= M(M* (k,z.). 
S(k,z) = M(k,z)<f) (k,z*), 



g*(z)$ (k, z») 



(25) 
(26) 



where g*(z) oc (1 + z)D(z) is the potential growth factor 
normalized to unity at last scattering (at which in our 
convention /nl is defined). 



In the peak-background split, we first divide any per- 
turbation .A(x) into a short- and long-wavelength piece: 



A(x) = A s (x)+A ; (x). 



(27) 



Here, the "long" wavelengths are those comparable to 
the distances over which we measure correlations. The 
"short" wavelengths are on much smaller scales, of or- 
der the Lagrangian extent of halos and less, whose peaks 
determine the locations of halos. For the Gaussian po- 
tential 0o, the two pieces <f>i = 4>o,l and 4> s = 0o,s are 
uncorrelated. The non-Gaussianity of Eq. (2) however 
induces a coupling of long- and short- wavelength pertur- 
bations in the potential as well as the matter density. 
Using Eq. (1) and the Poisson equation, Eq. (24), we can 
write the density <5 in terms of the Gaussian long- and 
short- wavelength pieces, 

<5(x) = ^(x) + <5 s (x)+2/ NL J d 3 x J d 3 yW (y,z) 

x j<Mx + y) [Si(x + z)+S s (x + z)] 

+</> s (x + y) [^(x + z) +(5 s (x + z)] 
1. 



- V x cj)o(x + y) • V x 4>o(x + z) 



(28) 



We are interested in the density field at the positions of 
density peaks. In those locations, we expect terms of the 
form (j>o 5 to be much larger than (V0o) 2 j since 5 is at 
least order unity, and additionally V0o is suppressed if 
we associate peaks in the density field with peaks in the 
potential (where V0o = 0). 

Next, we are mainly interested in the non-Gaussian 
effects on small-scale modes S s , since those result in the 
leading effect on halo formation. Separating the purely 
long- wavelength parts, we have 



6 a (x) = 5 s (x) + 2/ NL / d 3 x / d 3 yW (y,z) 



i(x + y)<5 s (x + z) (29) 

+0 s (x + y) [<y,(x + z) + 5 s (x + z)] 

We will now neglect the terms in the last line: 4> s 5i is very 
small, since 4> s ~ </>;, and 5i <C (J s ; and s <5s, while large, 
adds a purely stochastic contribution on large scales. We 
then obtain the final expression for the effect of primor- 
dial non-Gaussianity on small-scale modes in the context 
of the peak-background split: 

<5 s (x) =<5 s (x) (30) 
+2/nl f d 3 x J d 3 y Wo(y,z)<Mx + y)<5 s (x + z). 
Equivalently, in Fourier space: 

(31) 



5 s (k s ) = 6 s (k s ) 
d 3 k 



-2/] 



W (k,k s -k)<Mk)5 s (k s -k). 



5 



Note that by assumption, k s 3> k: the effect on the small- 
scale modes is determined by the behavior of the bispec- 
trum in the squeezed limit [Eq. (6)]. In the following, 
we denote short- wavelength modes with k s , while k will 
stand for long-wavelength modes. We reiterate that even 
for the local model, the coupling between density and 
potential is non-local, since Wq ^ 1. 



IV. FIRST-ORDER, NON-LOCAL HALO BIAS 

We assume that the number density of halos, per unit 
logarithmic mass, smoothed over a region of size A > fcj 1 
centered at position x is a function, 



ditions. It evaluates to 



b? 



(M, z) 



1 dn(M,z;p, P) 
n 91np 



1. 



(34) 



The second term in the last line of Eq. (33) arises if 
there are non-Gaussian initial conditions that correlate 
long- and short-wavelength modes. In this case, a change 
in <5;(k) induces changes in P(k s ) for each k s (hence the 
sum in Eq. (33)) which then induce changes in n(k). 

To evaluate the partial derivatives dP(k s )/dSi(k), we 
multiply Eq. (30) by <5 s (x'), letting r = x x', and take 
the expectation value. The first term that arises is the 
usual Gaussian two-point correlation function £ s (r). The 
non-Gaussian correction — that which correlates short- 
and long-wavelength modes — is then 



n^M,z;p[l + 6i(^)],P(k s ,Si(x))J, (32) A £,(r) = 4 | NL J d s y J £ z w ^ y> z )^( y )^ fl ( r + z ), (35) 



where M is the halo mass, z the redshift, p [1 + <5/(x)] 
is the density of matter in that region, and P(k s , <5;(x)) 
is the small-wavelength (linear) matter power spectrum 
in that region. In a smooth Universe — i.e., one with no 
long- wavelength fluctuations, 5i = 0, the halo abundance 
is the same everywhere. But if 5i ^ 0, then there will 
be fluctuations in the halo abundance with long-range 
correlations determined by long-range correlations in the 
matter density. Note that if primordial perturbations are 
Gaussian, then P(k s ) is the same everywhere. If, how- 
ever, primordial perturbations are non-Gaussian, and if 
that non-Gaussianity couples short- and long-wavelength 
modes of the density field, then the small-scale power 
spectrum P(k s ,5i(x)) may vary from one point in the 
Universe to another, as written in Eq. (32). 

Biasing describes the relative clustering of halos and 
matter. The Lagrangian bias is defined to be the 
ratio of the fractional halo- and matter-density pertur- 
bations. More precisely, here we will consider the scale- 
dependent bias 6i(fc) which describes the relative am- 
plitudes of Fourier modes of the halo abundance and 
matter density. To proceed, consider a single Fourier 
mode <5;(x) = (5;(fc)e lk x of the density field with wavevec- 
tor k and amplitude <5/(k). This Fourier mode in- 
duces a variation Sn(x) = <5n(k)e ik ' x with amplitude 
5h(k) = (dn/d(5(k)) x= o,(5,=o 5i(k). The scale-dependent 
Lagrangian bias is thus, 



b L {M,z;k) 



5h(k)/n dlnn(k) 



Sp/p 
ainn(k) 

a*i(k) 



dSiQs.) 

dlnn(k) dP(k s ) 

dp{k s ) a<5 ; (k) 



(33) 



E 



Here, the wavenumbers k of interest are for long- 
wavelength modes, those over which halo clustering is 
measured, and we have used the chain rule in the last 
line. 

The first term in the last line of Eq. (33) is the stan- 
dard result, that obtained assuming Gaussian initial con- 



where we have set x = without loss of generality. We 
now use 



A&(r) 



(27T) 



AP(k s ). 



(36) 



Expressing all quantities by their Fourier transforms, and 
performing the y, z integrals yields, 



A { , W=4AL /|^.</igL 



x WoQl^K - ki)M*i)P{k a )- 



(37) 



Hence, via Eq. (36), the effect of any quadratic primordial 
non-Gaussianity on the statistics of the short-wavelength 
modes is given by 



AP(k s ) = 4/, 



NL 



(2tt) 3 



W (k 1 ,k s -k 1 )^(k 1 )P(fc s 



(38) 



From this, we obtain the required partial derivatives, 

(39) 



9P(ks)_ M - l(k fAP(k s ) 



dS(k) 
dAP(k s 
dfa (k) 



4/ N LP(M^(k,k s -k), 



recalling that 4>i refers to the processed potential. Re- 
placing the sum ^ k by an integral J d 3 k s /(2Tr) 3 , we 
then obtain for the second term in the last line of Eq. (33) 
the non-Gaussian correction to the Lagrangian bias, 

Ab L (M,z;k) = 2/ NI >T 1 (fc) 

d 3 k* d Inn 



(2tt) 3 3P{k s ) 



Wo(k, k s -k)P(fc s 



(40) 



Thus, in general, non-Gaussianity changes the shape 
of P(k s ) as well as the amplitude. However, Eq. (40) 
simplifies for a very broad class of mass-function models, 
which only depend on the power spectrum through an 
overall normalization. We will turn to this case next. 
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A. Universal mass functions 

Most of the commonly used parametrizations of the 
halo abundance n(M, z\ ~p, P) are so-called universal mass 
functions [11, 29]. They can be written as 



(41) 



and thus only depend on cosmology through the back- 
ground density of matter p and the quantity, 



(42) 



where S c w 1.686 is the linearly extrapolated collapse 
threshold. Here, u\ is the variance of the density field 
on a scale R related to M by M = 47r/3pi? 3 , 



d 3 k 
J2nf 



F 2 R (k)P(k), 



(43) 



where Fr is the Fourier transform of a real-space tophat 
filter. The parametrization, Eq. (41), is sufficient to ob- 
tain the Gaussian and non-Gaussian halo bias. 

First, for the Gaussian halo bias we use Eq. (34). The 
effect of a long- wavelength perturbation Si on the collapse 
threhold is S c — > S c — Si , as derived in Refs. [30, 31]. This 
gives (see also Ref. [17]), 

b ^ = dlnn(p, v) = 2v d\i\{vf{v)) 



d5i 



dv 



On the other hand, it is clear from Eqs. (40) and (41) 
that the non-Gaussian halo bias will only enter through 

= _ 2 /Hvf{y)) = hS (45) 

a In (7r dv 
Next, we again use the chain rule, and, from Eq. (43), 
da 2 



dP(k s ) 



(46) 



We then have 



Ab L {k) = 2/ NI >T 1 (fc) 7 ^^#, (47) 
d 3 fc. 



cr w (k) 



(2tt) : 



-P(k s )F R {k s )W (K k s -k) 



The prefactor M^ 1 scales as k~ 2 . However, in general 
cr^, which involves Wo(k s ,k), will depend on k as well. 
Making use of the bias relations derived above, we have 



Ab L (k) = 2/ N i>r 1 (fc)M, 



(48) 



where, to clarify, the matter power spectrum P(k) is re- 
lated to the processed-potential power spectrum P$ (k) 
through P{k) = [M{k)] 2 P^ (k). 
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FIG. 1: Squeezed limit of W(q,[i) as function of q for equi- 
lateral (top) and folded (bottom) bispectra. Here, we have 
taken out the leading q scaling and adopted a fixed value of 
fi — 0.5. Shown is the exact result from Eq. (6), and the se- 
ries expansion of the squeezed limit to 0(q 2 , e 2 ) from § V. We 
also show the scale-free limit of this series expansion, which 
in general is not a good approximation. 



Let us first look at the local case. We use the squeezed 
limit [Eq. (23)] and W loc = 1. In this case, the change in 
the small-scale power spectrum is independent of fc s , so 
that only the overall amplitude is affected. We obtain: 



A^ ocal (fc) = 2/ NL 



M{k)T(k) 



(49) 



exactly matching the results of Refs. [17]. We will quan- 
tify the regime of validity of the squeezed limit in § VI. 

For a given bispectrum, Eq. (48) can be evaluated ex- 
actly in a straightforward way. However, it is also in- 
teresting to look at the scaling of the bispectrum in the 
squeezed limit, which allows for a simpler estimate of Ab. 



V. SQUEEZED LIMIT AND 
SCALE-DEPENDENCE OF BIAS 

As we have seen, the effect of primordial non- 
Gaussianity on the halo bias depends on Wo (k, k s — k) , 
where k s 3> k. This is known as the squeezed limit, 
since it corresponds to triangles where two sides are much 
longer than the third. In this Section, we derive the 
primordial kernel H^(k, k s — k) in the squeezed limit to 
show analytically how the non-Gaussian bias will behave 
for the various bispectra we consider. To leading order, 
the processed kernel is then given by Eq. (23). For our 
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numerical results in the next section, wc will use the 
full kernels, Eqs. (21)-(22). For definiteness, we assume 
a flat ACDM cosmology throughout, with n s — 0.958, 
n m = 0.28, h = 0.72, and er 8 = 0.8. 

Let us define q = k/k s and fx = k • k s . We then write 
the kernel as a power series in q. For example, in the 
simplest case, the local model, the squeezed limit is just 
W = 1. For the other bispectra, we use the k scaling of 
P*, 



P$(k) oc k* 



,-3+e 



(50) 



where e ~ —0.04 is the deviation from scale-free initial 
perturbations. The squeezed limit of the non-local bis- 
pectrum shapes introduced in § II can be easily derived in 
the limit of scale-free initial conditions. However, we will 
see that it is important to take into account a non-zero 
n a — 1. Since the prefactor in Eq. (40) goes as fc~ 2 oc q~ 2 , 



we will expand W up to terms oc q 2 
to second order in e. 



We will also expand 
For the equilateral model, this gives, 



W"*(q,n) 8qu H zed 3 g 2 ( 1 - M 2 



2s 



(lng-2 M 2 ) 



+ ¥ (21n 2 q- M 2 )) 



scale-free Q 2 /1 2\ 

3q ). 



(52) 
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FIG. 2: Non-Gaussian correction to the halo bias Ab^ calcu- 
lated in the peak-background-split (PBS) and local-bias ap- 
proaches as a function of k, for the equilateral (top) and SOSF 
(bottom) bispectra. We have assumed /nl = 100 in both 
cases and M — 10 13 M^/h halos at z = 0. We also show the 
scale- free approximations Eqs. (58)-(60). 



The squeezed limit for the orthogonal model is 



3q l + |(6 M 2 + /i-6) 



Thus, in the scale-free case, W eql scales exactly as k 2 /k 2 
in the squeezed limit, canceling the k~ 2 prefactor in 
Eq. (40) and leading to a scale-independent correction to 

the bias. This is not quite true in reality, since Eq. (51) +-( — IngH — [— 12/i 2 — fx+ (12 — /i) h\q\) 

contains logarithmic corrections, which can be quite im- 
portant for q -C 1 (Fig. 1). Note also that W = for 
fj, = 1 in the scale-free case, but is only suppressed by 
2e(\nq - 2)/3 -0.18 for q = 0.01. 

For the folded model, we have, 



3 V ' 2 

e 2 /, 9 q 



In 2 q + | [12^ 2 + 2 M In q + (fi - 24) In 2 q] 

cale— free , n 2 ^2/ i/^2 



3q + 9q*--q i (ji + W). 



(55) 



(56) 



^^fr>l / \ squeezed 3 / / , 0, ,\ Q 

W lo \q, n) = -q 1 + (4/i 2 + /i - 4) | 



£(l ng+ |[8 / x 2 + /x + (//-8)lng]) 



— (in 2 q + | [8^ 2 + 2/tlng + (/i - 16) In 2 q] 



scale— free 3 



(53) 



(54) 



In the scale-free case, the folded model shows a scaling 
proportional to q, with corrections going as q 2 . Fig. 1 
shows that the scale-free case matches the true squeezed 
limit somewhat more accurately than in the equilateral 
case. 



This bispectrum will therefore give rise to a scale- 
dependent bias similar to that for the folded model, but 
roughly twice as large in magnitude and opposite in sign. 

For the SOSF model, the kernel reduces in the squeezed 
limit to 



W 



SOSF 



squeezed 



(q, fi) aH "= 56.26 q l ( 1 - - In q l-^\nq 



(57) 



The q 2 scaling implies that this model will behave simi- 
larly to the equilateral model, although the prefactor of 
56 indicates that the amplitude of the effect will be much 
larger for a given /nl- 

As a rough approximation to the exact result for the 
non-Gaussian A&, we can take the scale-free squeezed 
limit derived above and in addition ignore any /i- 
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dependence. Prom Eq. (23) and Eq. (48), we then obtain 

6/nl 



AC 1 



A6 SOSF 



Ab\ o1 



A6 orth 



(k) 



(k) 



(k) 



(k) 



M(k)T(k) 

112/nl 
M(k)T(k) 

3/nl 
M{k)T{k) 

-6/nl 
M{k)T(k) 



b L S c k l 



b L S c k 2 



b L S c k- 



b L S c k- 



wherc we have denned 



'R.n 



d 3 k 



k n P(k) F R (k). 



(58) 



(59) 



(60) 



(61) 



(62) 



Since M. 1 oc k 2 , we see immediately that Afo cql and 
A5 sosf 

are expected to be roughly scale-invariant, while 
Ab and A6 orth are expected to scale as k^ 1 on large 
scales. Fig. 2 shows the exact PBS prediction (in the 
case of a universal mass function) for the non-Gaussian 
correction A6^ to the halo bias as a function of k, for 
the equilateral and SOSF models. We also show the ap- 
proximate relations assuming scale-free initial conditions 
derived in § V. We see that this approximation is only 
good to a factor of a few for these models. More impor- 
tantly, it docs not predict the remaining k dependence 
of the non-Gaussian correction, which is crucial to dis- 
entangle this effect from the ordinary Gaussian galaxy 
bias. The corresponding results for the folded and or- 
thogonal models are shown in Fig. 3. Here, the scale-free 
assumption gives a prediction accurate to ~ 20% on large 
scales. 



VI. PEAK-BACKGROUND SPLIT VS. LOCAL 
BIASING 

The effects on the halo power spectrum from non- 
local non-Gaussianity have also been derived in a local- 
Lagrangian-biasing (LLB) scheme [22, 23]. In this ap- 
proach, halos arc identified with high-density regions in 
the initial non-Gaussian density field. Specifically, the 
density <5r smoothed on a scale R (again related to the 
halo mass by M = 4ir/3pR 3 ) is required to be greater 
than the linear collapse threshold S c . This is equivalent 
to a local-biasing scheme relating the halo density field 
Sh to the smoothed (non- Gaussian) matter density field: 



4(x) = bJ R (x) + ^^(x) 



(63) 



where bz 2 is the second-order Lagrangian bias. In the 
high-peak limit which we will assume here, &£ = S c /a R , 
and bh,i = b\ (see discussion below). 

Since primordial non-Gaussianity induces a non-zero 
three-point function in the density field, the large-scale 
halo correlation function is modified to 
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FIG. 3: Same as Fig. 2, but for the folded (top) and orthog- 
onal (bottom) models. We have taken out a scaling factor of 
{k/O.Olh Mpc" 1 )- 1 . 



(3) 

where £r and £, R denote the two- and three-point func- 
tions of the smoothed density field, respectively. Note 
that, again, the three-point correlation function is eval- 
uated in the squeezed limit. In Fourier space, this reads 
[321 



P h (k)=blP R (k) + b L b L . 2 p S5S (k) 



p sss (k) 



d 3 k! 

(2^)3 



flfl(fc,fci,|ki-k|). 



(65) 
(66) 



Here, Pr(/s) = F R (k)P(k), and correspondingly for the 
bispectrum Br. Note that in this approach the correc- 
tion to the halo power spectrum comes about from the 
quadratic bias parameter. With this, we can write the 
effective correction to the linear halo bias as 



• SS5 (k) 



Ab ^ = l b ^ P -pW) 



1 p sss (k) 

b L6c 2 p 77 r . 67 



where for the second equality we have assumed the high- 
peak limit. We show the predictions of the local-bias 
model in Fig. 2 for the equilateral and SOSF models. We 
see that on the largest scales, the peak-background split 
prediction agrees precisely with that of local biasing. On 
the other hand, on intermediate scales k > 0.01 ft, Mpc -1 , 
the two approaches begin to diverge, with the peak- 
background split generally predicting lower halo bias cor- 
rections. We will now show how this comes about. 

We can write the bias correction in both the PBS and 
LLB schemes as (sec also [201) 



&(r) = bU R (r)+b L b L , 2 4 3) (x 1)Xl ,x 2 



(64) 



Ab L {k) = 2f Nh b L S c afT(k), 



(68) 
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with 

^PBS 



T 



LLB 



(A;) 



(k) 



d 3 h 



M{k)F R {k)J (2tt) 

1 r d 3 ki 
4hLP R (k) J Wf 



3 P fl (A;i)Wo(k,k 1 -k), 
(69) 



B«(A!,fci,|ki-k|; 



(70) 



The additional F^ 1 prefactor in J 71363 is necessary to 
match the local-biasing prediction, since the latter refers 
the halo bias to the smoothed density field [Eq. (63)]: do- 
ing the same in the PBS approach, via Eq. (33), Eq. (39) 
and using S R = F R 5, results in this prefactor. 

We now insert the definition of Wq [Eq. (21) or 
Eq. (22)] into Eq. (69), and use the relation between the 
smoothed density field and potential, 

P R {k) = M R {k)P* {k), 
B R (k!, fc 2 , fc 3 ) = M R (ki)M R (k 2 )M R (k 3 )B^ (ki,k2,k 3 ), 

where we have denned M R (k) = Ai(k)F R (k). It is then 
straightforward to show that 



^ Ub (k) = 

where 

A(k,kx) = 



1 



d 3 k! 

4/ NLj Pfl(fc)7 
x A(k,ki)C(k,ki), 



B R (k, fei, |ki - k|) 

(71) 



M R (h 



F R (h) 



M K (|ki-k|) Iki-kp^dki-kl) 1 



(72) 



and 



T 2 (|ki-k|) 

P*(|ki-k| 



x 1 + 



P*(|ki 



P*(fci) 



(73) 



(note that Pj> is the power spectrum for the primordial 
potential) if we use Wqi for the processed kernel, and 



C 2 (k,ki) 



2 1 



|ki 



P*o(*l) 



j^ (|ki-k|; 



(74) 



(and here P$ is the power spectrum for the processed 
potential) if we use W02 for the processed kernel. 

In the true squeezed limit, q = k/k\ becomes negligi- 
ble, so that ^4—5-1 and C\ — > 1 and C2 — > 1 (the sec- 
ond term in the brackets for both expressions becomes 
1, while the last term vanishes in this limit). Clearly, 
jrPBS _ jrioc - m £ u j g v] m j^ w hich explains why both ap- 
proaches agree on large scales. Going back to Eq. (69), we 
see that the integrand peaks where P R (ki) peaks, roughly 
around k\ ~ 0.02 h Mpc -1 . Furthermore, as k\ becomes 
comparable to k, |ki — k| can assume values much smaller 
than ki. Thus, we expect that deviations of J rPBS from 
jrLLB a pp ear a ^ least around k ~ 0.01 h Mpc -1 , which in- 
deed is visible in Fig. 2. Similar but somewhat smaller de- 
viations occur in the folded and orthogonal cases (Fig. 3), 
and even in the local model of non-Gaussianity (see 
Fig. 4). Hence, the agreement between the PBS and LLB 
approaches is restricted to the largest scales in the well- 
studied local model as well. However, since the halo-bias 
correction strongly declines towards smaller scales, the 
deviations between the two approaches have much less 
impact in the local model. We have confirmed numeri- 
cally that by introducing the correction factors A -1 and 
C -1 into the PBS expression Eq. (47), we recover the 
LLB prediction exactly. 

On a more physical level, the differences in scale- 
dependence between the PBS and LLB predictions come 
from the assumption of a separation of scales between 
the large-scale modes modulating the clustering of halos, 
and the small-scale modes that govern their formation. 
In the course of this assumption, several terms in the 
relation between short-wavelength and long-wavelength 
perturbations have been neglected (see § III). This as- 
sumption breaks down on smaller scales, however; modes 
with k ~ 0.01 h Mpc -1 , for example, contribute signif- 
icantly to the density variance o~ R , and if we want to 
calculate the clustering of halos at that k, the separation 
of scales does not hold anymore. On the other hand, no 
such assumption of separation of scales is made in the 
LLB approach. Here, the assumption is that halo forma- 
tion is purely a function of the local physical density. Of 
course, both approaches need to be benchmarked with 
simulations [33] in the end, which will have the final say 
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on the scaling of the halo bias in models with non-local 
non-Gaussianity. 

Even on large scales, there is a difference between 
the PBS and LLB predictions: in the PBS approach, 
Ab oc er^ dhxn/dlnaji, which for any universal mass 
function equals bi,&cla\- On the other hand, LLB pre- 
dicts that Ab oc 6l,2- While this agrees with the PBS 
prediction in the high-peak limit of the Prcss-Schechter 
theory, where bL.2 = b\ = b^Sc/a^, this is not the case in 
general. For example, the quadratic bias derived from the 
Sheth-Tormen mass function [31] in the high-peak limit is 
k>L,2 = Qb^c/cRi thus deviating from o~^d In nj d In an 
by a factor of q « 0.75. This number will be different 
in other mass function prescriptions. Also, the quadratic 
bias parameter has to change sign at a finite v, while 
6l<5 c /(t^ is always positive. Hence, outside of the high- 
peak limit the differences between LLB and PBS predic- 
tions can become larger. It is worth pointing out that 
since there is no prescription for small-scale physics con- 
tained in the local biasing approach, there is no analogous 
expression or prediction for the quantity dlnn/daR. In 
the end, a comparison with simulations will have to de- 
termine whether the non- Gaussian bias correction follows 
a scaling with dlnn/do-R, as predicted by the PBS, or 
with the quadratic bias bL,2, as predicted by local biasing 
- both of these predictions are falsifiable. 



VII. CONCLUSIONS 

The principal goal of the theory of large-scale struc- 
ture is to map the initial linear density field into the 
late-time non-linear matter density. One of the most 
fascinating recent applications of this theory is the ef- 
fect of primordial non-Gaussianity on the clustering of 
dark matter halos. Here, we derived predictions for gen- 
eral quadratic non-Gaussianity in the peak-background- 
split (PBS) approach, focusing on the case of universal 
halo mass functions. This complements the derivation 
of Ref. [22] in the local-bias model. We also show how 
the non-Gaussian effect on halo bias can be calculated 
beyond universal mass functions, by taking into account 
all dependences of the halo abundance on the statistics 
of small-scale fluctuations [Eq. (40)]. 

While we have only considered general quadratic non- 
Gaussianity here, it is possible to generalize these re- 
sults to the case of cubic non-Gaussianity, described by 
the trispectrum rather than the bispectrum of primordial 
fluctuations. For example, for the local cubic model, we 
expect the non-Gaussian halo bias to scale as ~ k~ 2 like 
in the local quadratic model [34]. A detailed investiga- 
tion of cubic non-Gaussianity is beyond the scope of this 
paper however. 

We show that both PBS and local bias approaches 
agree in their predictions on large scales, k < 
0.01 h Mpc -1 , when the high-peak limit is assumed in the 
local bias approach, but predictions diverge on smaller 
scales. The disagreement comes about from the break- 



down of the separation-of-scales assumption inherent in 
the PBS approach. Even on large scales, however, both 
approaches do not in general agree: in the PBS approach 
(with a universal mass function), the non-Gaussian effect 
on the halo bias scales with d\xin/d<jR = br,S c /<7ii, which 
quantifies the dependence of the halo number on the vari- 
ance of the small-scale fluctuations. By the nature of the 
approach, there is no analogous quantity in the local bias 
framework. Here, the halo bias correction scales with the 
quadratic bias parameter b^,2 which in general is not di- 
rectly related to dlnn/daR. This disagreement will also 
have consequences for predictions of higher order statis- 
tics (e.g., the bispectrum) of halos, where both bh,2 and 
the non-Gaussian bias contribution appear. 

N-body simulations will have the final say on which 
of the two approaches has the right scaling with halo 
mass (or equivalcntly, significance v), and with scale; it 
is possible that PBS fares better in the former respect, 
while LLB is closer to reality in the latter. 
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Appendix A: The local-model bispectrum from a 
nonlocal model 

In § II, we showed that the bispectrum _B,j>(fci, &2, fe) 
does not uniquely determine the kernel W(ki,k2). Sev- 
eral different kernels may give rise to the same bispec- 
trum. We moreover constructed explicitly two different 
processed kernels Wo(ki,k2) for a given bispectrum. Let 
us now use Eq. (22) for the processed kernel for the (pro- 
cessed) bispectrum for the local model. If we then invert 
Eq. (17) to obtain the primordial kernel that corresponds 
to this processed kernel, we find, 



W(ki,k 2 ) 



[TjkjTfa)] 2 (P1P2 + 2 perms) 
[T(fc 1 )T(fc 2 )] 2 J Pi J P2+2perms ' 



(Al) 



Upon substituting into Eq. (5), we recover the local- 
model bispectrum, Eq. (7). 

Moreover, the function T{k) that appears in Eq. (Al) 
need not be the transfer function for the local-model bis- 
pectrum to be recovered. Indeed, any kernel, 

m )= mm^lp^±i^, (A2) 

[/(fci)/(fc2)] 2 J PiP 2 + 2pcrms 
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for any function f(k) will yield the local- model bispcc- 
trum. And for any function f(k) that is not simply a con- 
stant, the configuration-space kernel W(y, z) obtained 
through the inverse transform of Eq. (3) will be nonlo- 
cal. _ 

In fact, the ambiguity in the kernel W is much more 
general: for any positive function of three arguments, 
g(ki, k 2 , k 3 ), and any bispectrum _B$, we can construct a 
kernel satisfying Eq. (5): 



W(k u k 2 ,k 3 ) 



g(ki,k 2 ,k 3 ) 



g(k\,k 2 , k 3 )P\P 2 + 2 perms 



2./'] 



NL 



(A3) 



We have therefore shown, by explicit construction, that 
there are models in which the potential $ is a nonlo- 
cal quadratic function of a Gaussian field <f> that have 
the same bispectrum as that for the local model. We 
therefore conclude that measurement of the local-model 
bispectrum does not necessarily imply that the potential 
has the local-model form given in Eq. (10). Note, though, 
that these different models may still differ in their pre- 
dicted C(/nl) trispectra. For now, we treat this simply 
as a mathematical curiosity, and we leave the investiga- 
tion of the implications of this result for model building 
for future work. 
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